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We present an analytical framework to study the escape rate from a metastable state under the 
influence of two external multiplicative cross-correlated noise processes. Starting from a phenomeno- 
logical stationary Langevin description with multiplicative noise processes, we have investigated the 
Kramers' theory for activated rate processes in a nonequilibrium open system (one-dimensional in 
nature) driven by two external cross-correlated noise processes which are Gaussian, stationary and 
delta correlated. Based on the Fokker-Planck description in phase space, we then derive the escape 
rate from a metastable state in the moderate to large friction limit to study the effect of degree 
of correlation on the same. By employing numerical simulation in the presence of external cross- 
correlated additive and multiplicative noises we check the validity of our analytical formalism for 
constant dissipation, which shows a satisfactory agreement between both the approaches for the 
specific choice of noise processes. It is evident both from analytical development and the corre- 
sponding numerical simulation that the enhancement of rate is possible by increasing the degree of 
correlation of the external fluctuations. 

PACS numbers: 05.40.-a, 02.50.Ey, 82.20.Uv 



I. INTRODUCTION 

Rate theory deals with the passage of a system from 
one stable local minima of the energy landscape to an- 
other via a potential energy barrier providing relevant in- 
formation on the long-time behavior of the system with 
different metastable states and hence is a very useful 
and pedagogical model for the microscopic description 
and understanding of a wide range of physical, chemical 
and biological phenomena. Examples include diffusion 
of atoms in solids or on surfaces, isomerization reactions 
in solution, electron transfer processes, and ligand bind- 
ing in proteins or protein folding reactions. Not only 
that the phenomena of flux transitions in superconduct- 
ing quantum interference devices can also be modeled by 
rate theory. Kramersi pointed out that the passage of 
a system from a metastable state due to thermal agi- 
tation can be considered as a model Brownian particle 
trapped in a one-dimensional well representing the reac- 
tant state which is separated by a barrier of finite height 
from a deeper well signifying the product state. The par- 
ticle is supposed to be immersed in a medium such that 
the medium exerts damping force on the particle but at 
the same time thermally activates it so that the particle 
may gain enough energy from its surroundings to cross 
the barrier. Thus, the surrounding fluid provides fric- 
tion and random noise. Conventionally, the friction is 
assumed to be Markovian and the noise, Gaussian and 
white. Kramers theory provides a useful framework for 
the explanation and prediction of rates and mechanisms 
for various barrier crossing phenomena and for chemical 
reactions in particular £iiL££JLLIL&12. 

Most of the above mentioned works are basically dealt 
either with phenomenological Langevin equation or by 
microscopic system-reservoir formulation with both bi- 



linear or non-linear coupling scheme. Although, the ef- 
fect of system-bath linear coupling and associate conse- 
quences in different branches of molecular sciences is well 
studied, the modeling and the proper microscopic inter- 
pretation of the nature of non-linear coupling and cor- 
responding manifestation on the various molecular phe- 
nomena is a very challenging task posed in front of re- 
searchers in this field. Tanimura and co-workers^ ex- 
plained the phenomena of elastic and inelastic relaxation 
mechanism and their cross effect in vibrational and Ra- 
man spectroscopy using non-linear system-bath model. 

It is now a well documented fact that solvent plays a 
very crucial role in computing the dynamics of a Brow- 
nian pAiticle^J^J^LlS^o The dynamics of a 
model Brownian particle in an inhomogeneous solvent 
generally leads to state dependent diffusion due to the 
appearance of multiplicative noise and state dependent 
dissipation! 19 i 20 i 21 Although modeling of Brownian mo- 
tion in presence of inhomogeneous solvent is a nontrivial 
task, several approaches have been made to study dif- 
ferent aspects of molecular sciences, e.g., activated rate 
processeS ) 5 i 21 ' 22 i 23 noise induced transport ) 19 i 20 i 24 ' 25 i 26 
stochastic resonance 2 ^ and laser and optics, 2 - 

In the overwhelming majority of the above men- 
tioned treatments, fluctuations experienced by the sys- 
tem of interest is of internal origin so that the fluc- 
tuations and dissipation get related through the cel- 
ebrated fluctuation-dissipation relation (FDR)* 2 . 9 How- 
ever, in some situations, the physical origin of dis- 
sipation and of fluctuations are different as well as 
independent i2£0 i 2L2&2L2Z&21i2Z. Thus, there exists no 
balance between the influx (efflux) of energy to (from) the 
system and hence there exist no FDR. As a result of this, 
such systems (usually termed as nonequilibrium open 
system 3 ^) do not approach thermal equilibrium asymp- 
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totically and many interesting phenomena are observed 
thereof-^ .^ 

To the best of our knowledge, creation of an open 
system in a nonequilibrium process can be realized in 
many ways. To mention, we briefly discuss about few of 
them in the following which are relevant to the present 
work. An additional perturbation applied to the heat 
bath can excite few bath modes which effectively leads 
to the creation of a nonstationary process where the 
nonstationarity mainly gets reflected in the dissipation 
kerneh 13 ' 14 ' 17 ' 37 Dynamics in this case is mainly guided 
by an irreversible process i 14 ' 17 An extensive analysis of 
this approach can be found in the very recent work of 
Popov and Hernandez^ 7 - External perturbation on the 
other hand can directly excite the system mode itself 
(leaving the bath mode alone) or the local bath modes 
(leaving the system mode alone). In both the situations 
the additional energy input by an independent source 
leads to a shift in the temperature, thus creating an ef- 
fective temperature like quantit y 19 ' 20 ! 21 ' 32 ! 33 ' 34 ! 38 keep- 
ing the underlying dynamical process stationary. An- 
other immediate effect of the breakdown of FDR in this 
case is the creation of a steady state instead of an equi- 
librium state in the long time limit. To study the escape 
rate from a metastable state under the influence of ex- 
ternal cross-correlated noise processes we have adopted 
one of the situations mentioned in the latter case to cre- 
ate an open system. In the present work we drive the 
reaction coordinate (the system mode) by multiplicative 
cross-correlated noise processes, leaving the bath modes 
as it is. In such a situation, the Brownian particle is en- 
ergized by an extra input of energy, in addition to the 
thermal energy provided by the heat bath. This leads 
to, what we have tried to study in this paper, an en- 
hancement of barrier crossing event in the activated rate 
processes within the framework of Markovian stationary 
dynamics. 

The barrier crossing dynamics with multiplicative and 
additive white noise processes aroused strong interest in 
the early eighties, where noise forces that are present 
simultaneously in the dynamical system were usually 
treated as random variables uncorrelated with each other. 
However, there are situations where fluctuations in some 
stochastic process may have common origin. If this hap- 
pens then the statistical properties of the fluctuations 
should not be very much different and can be correlated 
to each other. Nonlinear stochastic systems including 
noise terms have drawn interest on a wide scale owing 
to their numerous applications in the field of molecu- 
lar sciences. Generally it is observed that under such a 
situation, the noise affects the dynamics through a sys- 
tem variable. The cross-correlated noise processes were 
first considered by Fedchenia^ 9 - in the context of hydrody- 
namics of vortex flow where the author introduced cross- 
correlation among the fluctuations from a common origin 
that appear in the time evolution equation of dimension- 
less modes of flow rates. The interference of additive 
and multiplicative white noise processes in the kinetics 



of the bistable systems was first considered by Fulinski 
and Telejko^ - where they mentioned the physical possi- 
bility of a cross-correlated noise. Madureira et alM- have 
pointed out the probability of cross-correlated noise in a 
realistic model (ballast resistor) showing bi-stable behav- 
ior of the system. Recently, Mei et a/4 2 - have studied the 
effects of correlations between additive and multiplica- 
tive noise on relaxation time in a bi-stable system driven 
by cross-correlated noise. 

It is now well accepted that the effect of correlation 
between additive and multiplicative noise (which is also 
termed as correlated noise processes) is considered indis- 
pensable in explaining phenomena such as steady state 
properties of a single mode laser^ 3 - bistable kinetics^ 4 - 
stochastic resonance in linear system,— steady-state en- 
tropy production^ stochastic resonance^ 7 - and transport 
of particles^ 8 - etc. ZhvM- investigated theoretically the 
statistical fluctuations of a single mode laser that in- 
clude correlations between additive and multiplicative 
white noises and showed that the effect of correlation 
can lead to larger intensity fluctuations. One can uti- 
lize this intensity fluctuations to induce a narrow peak 
in the probe absorption spectrum as well as significantly 
modify the emission spectra of matter strongly resonant 
with laser field. As shown by Ai et a/4 9 - one can use the 
logistic differential equation to analyze the effects of en- 
vironmental fluctuations on cancer cell growth using cor- 
related Gaussian white noise scheme whereby the degree 
of correlation of the noise (environmental intensive fluc- 
tuations) can cause tumor cell extinction. Berdichevsky 
and Gittermai>S showed that the maxima of signal to 
noise ratio, as a function of the asymmetry of noise, dis- 
appears in the absence of the coupling between additive 
and multiplicative white noises. Very recently Ghosh et 
al^ have used multiplicative correlated noise processes 
to study the splitting of Kramers' rate in a symmetric 
triple well potential. 

As the presence of the cross-correlated noise changes 
the Langevin dynamics of the system ; 38 ' 51 it is expected 
that there may exist some additional effect of cross- 
correlation on the escape rate of a metastable state. 
Study of this additional effect can perfectly serve as a 
motivation of our work presented in this paper. To 
achieve this we extend our recently developed theoreti- 
cal approach^ 8 , to study the effect of multiplicative cross- 
correlated noise on the escape rate from a metastable 
state. In this present study, the external fluctuation ap- 
plied to the system under consideration is assumed to 
be independent of the system's characteristic dissipation. 
In this article, we study the reaction rate (in one di- 
mension) under the influence of both the internal (ther- 
mal) noise and external (non-thermal) cross-correlated 
noises simultaneously to examine the role of correlation 
between external delta correlated fluctuations on the es- 
cape rate from a metastable state when the dissipation is 
space dependent and the noises appear multiplicatively 
in the dynamical equation. In our model the system 
is externally driven by two cross-correlated fluctuations. 



3 



Starting from phenomenological Langevin equation with 
space dependent dissipation and multiplicative noises, we 
construct the corresponding Fokker-Planck equation in 
phase space. Under proper boundary conditions, wc solve 
the Fokker-Planck equation to calculate the escape rate 
for moderate to large dissipation. In the numerical im- 
plementation of our development the steady state rate 
expression (derived analytically) is checked with stochas- 
tic simulation to establish the fact that the incorporation 
of external perturbation through cross-correlated fluctu- 
ations in to the traditional Kramers' model enhances the 
rate across the barrier top. 

The organization of the paper is as follows. In section 
II, starting from a phenomenological Langevin equation 
for an open system nonlinearly coupled with the envi- 
ronment and simultaneously acted upon by two cross- 
correlated multiplicative white noise processes, we con- 
struct the Fokker-Planck description of the underlying 
stochastic process which is multiplicative in general. We 
then calculate the escape rate from a metastable state to 
examine the barrier crossing dynamics. A typical exam- 
ple has been considered in section III to study the effect of 
the cross-correlated fluctuations on the escape rate. The 
paper is then concluded in section IV. To make the paper 
self contained we provide the derivation of the Fokker- 
Planck equation and the calculation of escape rate in the 
two appendices. 



Ub is the Boltzmann constant. e(t) and ir(t) are Gaussian 
white noise processes with statistical properties 



(e(t)) 
(e(t)e(t')) 

(e(i)Tr(i')) 



W*)> = o 
2D e S(t - t') 
2DJ{t - t') 
(n(t)e(t>)) = 



(3a) 
(3b) 
(3c) 

2X^D e D w 5(t - 1'). (3d) 



In the above equations (|3bM3dp . D e and are the 
strength of the fluctuations e(t) and 7r(i), respectively 
and A (0 ^ A < 1) denotes the degree of correlation 
between the noise processes e(t) and 7r(i). The internal 
noise f(t) is also assumed to be Gaussian and delta cor- 
related with statistical properties 



(/(*)) =0 and (f(t)f(t'))= 2S(t-t'). 



(4) 



In the above equations (••■) implies the average over 
the realizations of the noise (external or internal) pro- 
cesses. Eq.([2]) along with the second relation of Eq.((4]) 
comprises the fluctuation-dissipation which relates the 
damping function T(x) with the fluctuation f(t). The 
external noise processes e{t) and ir{t) are independent of 
the dissipation function and there is no corresponding 
fluctuation-dissipation relation. We further assume that 
f(t) is independent of e(t) and ir(t) so that 



II. CORRELATED NOISE INDUCED ESCAPE 
FROM A METASTABLE STATE 

We consider the motion of a particle of unit mass mov- 
ing in a Kramers type potential V{x) such that it is acted 
upon by random force /(<) of internal origin, i.e., f(t) 
originates due to the coupling of the system with its en- 
vironment and hence is connected to the friction through 
the fluctuation-dissipation relation. Apart from the in- 
ternal noise, we assume that the system is acted upon 
by two external Gaussian noise processes, e(t) and 7r(i), 
both of which have a common origin and consequently 
are correlated. Thus, from the very mode of description 
of our model, it is evident that the system is open in 
nature as it is driven externally by two correlated fluctu- 
ations. The dynamics of the particle is governed by the 
Langevin equation 

x = -T(x)x - V'(x) + h(x)f(t) + gi(x)e(t) + g 2 (x)n(t), 

(1) 

with 



</(t)e(t')> = (/(t)Tr(f )) = 0. 



(5) 



h(x) = \J feTr(i), 



(2) 



where T(x) is the space dependent friction that arises due 
to non-linear system-reservoir coupling.— Here, gi(x) 
and g 2 (x) are two arbitrary smooth functions of x and 
their presence makes the external noise processes multi- 
plicative. T is the thermal equilibrium temperature and 



In the absence of the external noise processes, the 
system being closed, the fluctuation-dissipation relation 
eventually brings it to a stationary state and conse- 
quently one can examine the barrier dynamics of the 
system using standard methods applicable for a ther- 
modynamically closed system. The external fluctuations 
and their correlation modify the dynamics of activation 
in two ways. First, they influence the dynamics in the 
region around the barrier top so that the effective sta- 
tionary flux across it gets modified. Second, in the pres- 
ence of these fluctuations, the equilibrium distribution 
of the source well is disturbed so that one has to con- 
sider a new stationary distribution, if any, instead of the 
standard equilibrium Boltzmann distribution. This new 
stationary distribution must be a solution of the Fokker- 
Planck equation around the bottom of the source well 
region and serves as an appropriate boundary condition 
analogous to Kramers problem. 

Keeping this in mind we write the corresponding 
Fokker-Planck equation in the phase space [i.e. in (x, v) 
space of the system, where v = x = dx/dt], describing the 
dynamics of the Langevin equation [see Appendix-A 
for detailed derivation] 



dP 
~dt 



dP 
dx 
+T(x)P, 



[T(x)v + V'(x)\ 



dP 

dv 



A(x) 



d 2 p 

dv 2 



(6) 
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where 



equation ([9]), namely 



A{x) = k B TT{x)+g 2 {x) 

g(x) = {D e gf(x) + D e D 7T g 1 (x)g 2 (x) 

+D 7I gl(x)} 1/2 . 



(7) 



(8) 



The diffusion coefficient A(x) in the Fokker-Planck equa- 
tion ([6]) is an implicit function of the correlation param- 
eter A (see the expression of g(x)). Thus by increasing 
the value of A one can increase the value of the diffu- 
sion coefficient. This property gets directly reflected in 
the final rate expression and in the expression of effective 
temperature (see Eqs.(12) and (13)). In deriving Eq.©, 
we have made use of Eqs. (|3aM3dl) and (j4|) and the fact 
that f(t) is independent of e(t) and Tt(t). It should be 
noted that when the noise is purely internal and for linear 
system- reservoir coupling, Eq.© reduces to the Kramers 
equation^ 

Kramers' model for a chemical reaction consists of a 
particle undergoing Brownian motion whose coordinate 
x corresponds to the reaction coordinate and v = dx/dt 
the reaction rate. In addition to that in Kramers' original 
treatment,— the dynamics of the Brownian particle was 
governed by Markovian random process. Since the work 
of Kramers, a number of workers have extended Kramers 
model for non-Markovian case and for state dependent 
diffusion to derive the expression for the escape rate£ In 
order to allow ourselves a comparison with the Fokker- 
Planck equation of other forms, we note that though the 
underlying dynamics is Markovian, the diffusion coeffi- 
cient in Eq.® is coordinate dependent. It is customary 
to get rid of this dependence by approximating the co- 
efficients at the barrier top or potential well where we 
need the steady-state solution of Eq.©. One may also 
use mean field solution of Eq. ([8]) obtained by neglecting 
the fluctuation terms and putting approximate station- 
ary condition in the diffusion coefficient. 

For harmonic oscillator with frequency luqi V( x ) ~ 
Eq + uJq(x — xq) 2 /2, the linearized version of the Fokker- 
Planck equation can be represented as 



dP 
~dt 



dP 

-^ + r ( , o) P. 



[T(x )v 



■uj {x-x )} — 



+A 



d 2 P 
' dv 2 ' 



(9) 



where A = kBTT(x ) + g 2 (x ), is calculated at the bot- 



tom of the potential (x w xq). 
solution of Eq. © becomes 



The general steady state 



V 2 L)l(x - X f 



2Dn 



2Dn 



(10) 



with Dq = Aq/T(xq) and Z is the normalization con- 
stant. The solution (fTU|) can be verified by direct substi- 
tution in the steady state version of the Fokker-Planck 



dP st 

-v^+T(x )P^ + [T(x )v 



-A 



dx 

Q2pst 

dv 2 



-uj 2 (x 



o. 



It should be noted that the distribution p0[) is not an 
usual equilibrium distribution, rather it serves as sta- 
tionary distribution for the non-equilibrium open system. 
This stationary distribution plays the role of an equilib- 
rium distribution of the closed system which can, how- 
ever, be recovered in the absence of the external noise. 

We now embark on the problem of decay of a meta 
stable state in the presence of external cross-correlated 
noise processes. In Kramers' approach, the particle co- 
ordinate x corresponds to the reaction coordinate, and 
its values at the minima of the potential well V(x) are 
separated by a potential barrier to describe the reactant 
and product states. Linearizing the motion around the 
barrier top at i sj ij, the steady state version of the 
Fokker-Planck equation corresponding to Eq. ([6]) reads as 



fipst 

ur b 2 1 

-V— u b {x - 



-A 



dx 

Q2pst 



Xb) 



8P b st 
dv 



+ T(x b 



djvPf) 
dv 



dv 2 



0. 



(11) 



where, V(x) ~ Eb — u\(x — Xb) 2 /2 with lo 2 > 0, and the 
suffix '6' indicates that coefficients are to be calculated 
using the general definition of A at the barrier top. 

After imposing appropriate physically motivated 
boundary conditions we then derive the escape rate for 
nonequilibrium open systems valid in the moderate to 
strong friction limit [see Appcndix-B for detailed calcu- 
lation] 



A 



uj D b 



2tt i/D^ V 1 + AD. 



■ exp 



E 



(12) 



where E = E b — Eq is the potential barrier height, and 
Db = Ab/T(xb)- Already we have mentioned in the in- 
troduction, open system mean the system is not ther- 
modynamically closed as it is driven by two noises of 
external origin and consequently there is no fluctuation- 
dissipation relation which is mainly responsible to estab- 
lish the thermal equilibrium (for closed system). It is 
also important to note that the escape from a metastable 
state intrinsically a non-equilibrium phenomenon as fluc- 
tuations establish a steady mass motion to give rise to a 
non- vanishing constant current over the potential barrier. 
In the steady state rate expression (1121) . all information 
about the thermal temperature due to internal noise pro- 
cesses, the strength of both the external noises and the 
position dependent dissipation, T(x), are hidden in the 
quantities Do, Db and A. Furthermore, the rate constant 
k is an implicit function of the degree of correlation, A, 
between the external noise processes. We conclude this 
section by mentioning the fact that the A-dependence 
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of the effective temperature (£>&) of our model does not 
depend on the detailed forms of the coupling functions 
gi(x) and g 2 (x). 



III. ANALYSIS OF THE INTERFERENCE 
BETWEEN TWO EXTERNAL FLUCTUATIONS 

In this section we present the numerical implementa- 
tion of our presently developed theory to establish its 
potentiality and applicability to demonstrate the effect 
of correlation of two external noises e(t) and ir(t) on the 
rate. Eq. (fl~2")) is the theoretical expression for the escape 
rate and is applicable for any two functions g% (x) and 
g 2 (x) appeared in equation |T]). To test the validity of 
the rate expression (j 1 2[) . one needs to assume the par- 
ticular forms of gi(x) and g 2 (x). In addition to that 
an explicit form of the damping term Y(x) is needed. 
For simplicity we use a constant dissipation T(x) = 7 in 
our numerical simulation. Thus the explicit expression 
for the escape rate for constant dissipation and arbitrary 
gi{x) and g 2 (x) becomes 



2-KUJb 



x exp 



-fk B T + g 2 (x b ) 



lk B T 



jE 



1/2 



7' 



1/2 



ryk B T + g 2 (xb) 



(13) 



where the quantities g(xo) and g{xb) are evaluated 
around the potential minima (x ~ xq) and maxima 
(x ~ .Tfc), respectively, using the function g(x) defined 
in Eq.©. In the above equation the quantity g 2 (xb)/ r y 
in the exponential factor together with the thermal en- 
ergy k B T constitutes the effective temperature, a typical 
signature of the open system. 17 i 21 i 32 i 33 i 34 i 38 Another ef- 
fect of space dependent diffusion in the dynamics is the 
presence of reaction coordinate (xq and Xb) in the steady 
state rate expression (|13[) which is typically absent in 
the standard Kramers' expression 1 , (see Eq.(fT4"f below). 
For D e — and D, = 0, i.e., in the absence of exter- 
nal driving, Eq. (|13p yields Kramers' rate expression for 
pure thermal fluctuations valid in the moderate to strong 
damping limit 1 



ramers 



^0 

lltLOb 



x exp 



2 + 
E 



1/2 



(14) 



which can be also verified easily using the explicit forms 
of the parameters Do, Db and A in Eq. (fT2")) in the limit 
e(t) = Tr(t) = 0. It is thus clear that in the absence of the 
external noise processes, the derived rate expression (jTSJ) 
reduces to standard Kramers' escape rate in the moderate 
to large damping regime. 

To study the dynamics, we consider a model cubic po- 
tential of the form V(x) = b\x 2 — b 2 x 3 where, b\ and 




FIG. 1: Plot of barrier crossing rate k as a function of dissi- 
pation constant 7 for various values of correlation parameter 
A and for gi(x) = x and Q2(x) = 1. The solid lines are drawn 
from the theoretical expression, Eq. (|13p and the symbols are 
the results of numerical simulation of Eq.([T]). The values of 
the parameters used are fcsT = 0.1, D ( — = 0.1 and A 
0. 0.5 and 0.9. 



&2 are the two constant parameters with 61,62 > 0, so 
that the potential barrier height becomes 46 3 /27&2 and 
Xb = 26 2 /3&i. In our numerical simulation we have used 
b\ =62 = 1. While numerically solving the Langevin 
equation ([I]) for constant dissipation, 7, we used a specific 
combination of the multiplicative terms, e.g., <?i(x) = x 
and g2(x) = 1. Although both the external noise pro- 
cesses in our model are multiplicative in nature, the spe- 
cific choice of the external noise processes we adopt in 
the numerical simulation captures the essential feature 
of our model (see the discussion in the following para- 
graph). We then numerically solve the Langevin equa- 
tion {1} by employing stochastic simulation algorithm. 53 
The numerical rate has been defined as the inverse of the 
mean first passage tim o 21 ' 32 i 34 ' 54 and has been calculated 
by averaging over 10,000 trajectories. To ensure the sta- 
bility of our simulation, we have used a small integration 
time step At = 0.001. 

In Fig. 1, we show the comparison between the simu- 
lation result and theoretical expression where the escape 
rate k is plotted as a function of the dissipation constant, 
7, in the moderate to large damping regime where our 
theory is valid, for various values of the degree of corre- 
lation A and observe that for a given 7, the escape rate 
k increases with an increase in A. The result shows a 
reasonable agreement between the theory and the sim- 
ulation. It is easy to check that for a particular com- 
bination of g±(x) and gi (x), by increasing the correla- 
tion parameter A one basically increases the value of the 
function g(x) (see Eq.(|S])) evaluated at x i=s Xb- This in- 
crement in g{xb) increases the effective temperature like 
quantity k B T + (g 2 (xb) / '7) for a fixed value of 7 which in 
turn pumps more energy into the system which eventu- 
ally helps in crossing the barrier and increases the escape 
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rate. Another way to explain this phenomenon is to look 
at the diffusion coefficient, Af, across the barrier. From 
the expression of A(x) (see Eq.Q) evaluated at x « Xb, 
i.e., Af, — jIcbT + g 2 {xb) it is clear that g(xb) increases 
with an increase in A for a particular choice of 171 (x) and 
#2(2) and for a fixed value of "fknT. This increment in 
the diffusion across the barrier enhances the reaction rate 
as observed. This is the central result of this paper. 



IV. CONCLUSIONS 

In conclusion, we have extended our recently devel- 
oped theoretical approach of studying escape rate from 
a metastable state under the influence of external addi- 
tive cross-correlated noise processes^ to investigate the 
effect of multiplicative noise. In contrast to our previous 
approach of using effective correlated noise constructed 
from two additive colored noise processes (see Eq.(l) of 
Ref. 38), the correlated noise used in the present work has 
been constructed from multiplicative white noises. The 
multiplicative nature of the correlated noises introduces a 
space dependent diffusion in the resultant Fokkcr-Planck 
equation ^ and modifies the exponential as well as the 
pre-exponential factor of the steady state rate expression 
(fT3"]) . To check the validity of the steady state analyt- 
ical rate expression, we have performed numerical sim- 
ulation of the starting Langevin equation (TT]) with one 
multiplicative noise and the other additive one, which 
shows a satisfactory agreement between the theory and 
the numerics. The analytical development as well as the 
corresponding numerical simulation lead us to conclude 
that the enhancement of rate is possible by increasing the 
degree of correlation of the external fluctuations. So far, 
we have used the formalism of Markovian stochastic pro- 
cesses in this paper which can be extended to look into 
the dynamics within the framework of non-Markovian 
formalism. We plan to address this issue in our future 
communication. 
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APPENDIX A: DERIVATION OF THE 
FOKKER-PLANCK EQUATION 

In this appendix we give the detailed calculation of 
constructing the Fokker-Planck equation ^ with space 
dependent diffusion corresponding to equation (JTJ) which 



can be written as: 

x = v, 

v = -r(x)v-V'(x) + h(x)f(t) (Al) 
+ 9i(x)e(t) + g 2 (x)ir(t), 

We then rewrite the above equation in the following form: 

ui = F 1 (u 1 ,u 2 ,t; f(t),e(t),7r(t)), 

u 2 = F 2 (ui, «a, t; f(t), e(t), tt(*)), (A2) 

where we use the following abbreviations 

u\ = x, 1*2 = v, (A3) 
F 1 =v,F 2 = -r(x)v-V'(x)+h(x)f(t)+ 9l (x)e(t) 
+g 2 (x)Tr(t). (A4) 

The vector u with components u\ and u 2 thus represents 
a point in a two-dimensional 'phase space' and equation 
(|A2[) determines the velocity at each point in the phase 
space. The conservation of the phase points now asserts 
the following linear equations of motion for density p(u, t) 
in phase spaced 

d d 

-p(u, t) = -J2 g^ F n(u, t, /(*), e(t), 7T(t)M«> *)> 

n—1 n 

(A5) 

or more compactly 

| = -V-FP. (A6) 

Our next task is to find out a differential equation whose 
average solution is given by (/?),— where the stochas- 
tic averaging has to be performed over both the inter- 
nal noise process f(t) and the external noise processes 
e(t) and n(t). To this end we note that V • F can be 
partitioned into two parts: a constant part V • Fq and a 
fluctuating part V • F\ (t) containing these fluctuations. 
Thus, we write 

V-F(u,tJ(t),e(t),7r(t)) 

= V • F (u) + /3V • Fi(u, t, f(t), e(t), n(t)), 

where /3 is a parameter (we put it as an external param- 
eter to keep track of the perturbation equation; we put 
j3 = 1 at the end of the calculation) and also note that 
(Fi(t)} = 0. Equation (|A6|) therefore takes the following 
form: 

P (u,t) = (A +(3A 1 )p(u,t), (A7) 

where Ao = — V • Fq and A\ = V • F\. The symbol V is 
used for the operator that differentiates everything that 
comes after it with respect to u. Making use of one of the 
main results for the theory of linear equation of the form 
(IA7[) with multiplicative noise^ we derive an average 
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equation for p [ (p) — P(u,t) the probability density of 
u(t)], 

BP ( f°° 

= \A o +0* J d T (A 1 (t)eMrA )A 1 (t-T)) 



x exp(-rA )}P. 



(A8) 



Equation (|A8|) is exact when the correlation times of fluc- 
tuations tend to zero. Using the expressions for Aq and 
A\ we obtain 



dt 



pOO 

F + (3 2 dr{W-F 1 {t) 
Jo 



x exp(-rV • F )V ■ F x {t - r)) exp(rV • F )} P. 

(A9) 

The operator exp(rV-i 7 b) in the above equation provides 
the solution to the equation 



dy(u,t) 
dt 



= -V-F y(u,t), 



(A10) 



(y signifies the unperturbed part of p), which can be 
found explicitly in terms of characteristic curves. The 
equation 



u = F (u), 



(All) 



for fixed t determines a mapping from u{t = 0) to u(t) 
i.e., u — ► u T with the inverse (u T )~ T = u. The solution 
of equation (|A10|) is 



y(u,t) = y(u *,0) 



d(t»-*) 



d{u) 



exp(-fV-F )y(M,0), 



(A12) 

\d(u t )/d(u)\ being a Jacobian determinant. The effect 
of exp(— tV ■ Fq) on y(u) is as follows: 



exp(-tV.foM«,0) = v(«- t ,0) 



d(u 



d(u) 



(A13) 



The above simplification when put in equation (|A9 
yields 



dP{u,t) 
dt 



= V • <{ -F + [3 I dr 



d(u- T ) 



d(u) 

x(F 1 {u,t)V^ T -F 1 {u- T ,t-T)) 
d(u) 



d{u- T ) 



P(u,t). 



(A14) 



where V_ T denotes differentiation with respect to u~ T . 
We put (3 ~ 1 for the rest of the treatment. Identifying 
u\ = x and u 2 = v we now write 

F i = v : F u = 0, 

F 02 = -T(x)v-V(x), 

Fi 2 = h(x)f(t)+ gi {x)e(t)+g 2 (x)iT{t). (A15) 



In this situation, equation (|A14p now reduces to 



dP 
~dt 



+ fo J dT d h (x).f(t) + 9i(x)e(t) + g 2 (x)7r(t) 



+ -^{h(x- T )f(t -t)+ g x {x- T )e(t - t) 



+ 9 2(x~ T Mt-T))}])P, 



(A16) 



where we have used the fact that the Jacobian obeys the 
equation 



dt 



log 



d{x\v t ) 



d(x, v) 



| + « { _r. + v W} -r, 



so that the Jacobian becomes exp(— T(x)t). Now neglect- 
ing the terms 0(t 2 ) we may have x~ T = x — tv and 
v ~ T = v + Ttv + tV{x). The above two equations yield 



dv 



d _ ^ . d ^ ^_ d 
T dv dx 



(A17) 



Taking this in to consideration of equation (|A17j) . equa- 
tion (TAT6|) can be simplified in the following form: 



dP 
~dt 



d(vP) d 



dx 
d 2 



dv 



T[x)v + V'(x)]P 



+ — [k B TT(x)+g 2 (x)}P, 



(A18) 



where 



g{x) = [D e gf(x) + 2X^ D e D^ gi (x)g 2 (x) 



+D 7T g 2 (x)} 



1/2 



Defining 



A{x) = k B TT(x)+g 2 (x), 
the above equation (|A18[) can be written as 



(A19) 



(A20) 



+T(x)P, (A21) 
which is our required Fokker- Planck equation. 



APPENDIX B: CALCULATION OF THE ESCAPE 
RATE 

In this appendix we show the detailed calculation 
of the escape rate for nonequilibrium open system. 
The technique we adopt here resembles our previous 
approache o 32 i 33 i 38 but makes the current paper self con- 
tained. 
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Following Kramers we make ansatz that the non- 
equilibrium steady state probability P h st (x,v) generat- 
ing a non-vanishing diffusion current across the barrier 
is given by 



Pg t (x,v) = exp[- — - 



D h 



G(x,v), (Bl) 



where D b = A b /Y(x b ). Inserting Eq.([B"l]) in ([TT]1. we 
obtain the equation for G{x, v) using the steady state in 
the neighborhood of x b 

dG r 9/ \ , , ,0G . d 2 G „ 
- v— - \u 2 (x - x b ) + T{x)v}— + A b — T = 0. (B2) 



dx dv 
At this point we set 

y = v + a(x - Xb), 



dv 2 



(B3) 



where a is a constant to be determined. With the help 
of the transformation (|B3|) . Eq. (|B2[) reduces to 

,dG 
dy 



A "^T - I* - x b) + {F(xb) + a}v}^ = 0. (B4) 
dy 1 



Now let 



\w\{x - Xb) + {T(x b ) + a}v] 



(B5) 



where [i is another constant. By virtue of the relation 
(|B5l) . Eq. (|B4|) becomes 



dy 2 dy 



(B6) 



where A = fj,/A b with A b = k B TT(x b ) + g 2 {x b ). The 
constant /i and a must satisfy the following equations 
simultaneously: 



fia — uj 2 and — fj, = T(x b ) + a. 



(B7) 



This implies that the constant a must satisfy the follow- 
ing quadratic equation 

a 2 + T(x b )a — uj 2 = 0, 

which allows the solution for a as 



Thus, the general solution of (|B6j) is 



G(y)=G 2 J exp(-— jdz + G u 



(B8) 



(B9) 



where G\ and G 2 are two constants of integration. We 
look for a solution which vanishes for large x. This con- 
dition is satisfied if the integration in (|B9|) remains finite 
for \y\ — > +oo. This implies that A > so that only o_ 
becomes relevant. Then the requirement P b (x, v) — > for 
x — > +oo yields 



G\ — G A 



Thus we have 

G(y) = G 2 
and correspondingly 

P b st (x,v) = G 2 

x exp 



exp 



dz 



(Bll) 



7T 

2A jo 

v 2 V{x) 
2D~ b ~ 



y / Az 2 

exp I — ] dz 



Di, 



(B12) 



The current across the barrier associated with this steady 
state distribution is given by 



J 



vP b st {x w x b , v)dv, 



which may be evaluated using (|B12|) and the linearized 
version of V(x), namely V(x) « E b — uj 2 (x — x b ) 2 /2 as 

E b 



G 2 



2tt 



— D b exp 



(B13) 



To determine the remaining constant G 2 we note that 
as x — * — oo, the pre-exponential factor in (|B12[) reduces 
to G 2 \/2-k I A. We then obtain the reduced distribution 
function in x as 



P?(x - -oo) = 2ttG 2 ^ exp (-^ 



(B14) 



where we have used the definition for the reduced 
distribution2£ii&2& as 



P(x) 



P{x, v)dv. 



Similarly, we derive the reduced distribution in the left 
well around x ~ xq using Eq. (|10[) where the linearized 
potential is V(x) w Eq + to 2 {x — xq) 2 /2, 



P^(x) = -V27rA)exp 



u> 2 (x - X Q f 
2D 



(B15) 



with the normalization constant given by 1/Z = 
u>o/(2irDo) The comparison of the distribution (|B14[) and 
(|B15[) near x sa xq, gives, 



G- 2 = 



D b 27rV27rA, 



(B16) 



Hence, from (|B13|) . the normalized current or the barrier 
crossing rate k, for moderate to large friction regime is 
given by 



A 



E 



~ 2 7 rVA;Vl + AD b CXP l v D b 



(B17) 



(BIO) where E = E b — Eq is the potential barrier height. 
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